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Abstract 

We  propose  a  remote-sensing  technique  to  measure  temperature  in  the  lower  thermosphere  with  a  resonant  Raman  lidar. 
A  ground-based  pulsed  laser  operating  at  630.0304  (636.3776)nm  excites  3P2  (3Pi)  multiplet  level  of  the  ground  electronic 
state  of  atomic  oxygen  in  the  atmosphere  to  the  electronically  excited  *D2  state  and  the  back-scattered  photons  at  636.3776 
(630.0304)  nm,  while  the  atom  transitions  to  3Pi  (3P2),  are  detected.  Using  the  backscattering  Raman  cross  sections 
calculated  here  we  show:  (1)  For  the  range  of  altitudes  in  the  lower  thermosphere  where  the  fine-structure  multiplets  of 
atomic  oxygen  are  in  thermodynamic  equilibrium  with  the  local  translational  temperature  (LTE)  and  the  electronically 
excited  intermediate  state  ‘D2  remains  relaxed  primarily  by  collisions  with  N2  and  02,  the  ratio  of  the  backscattered  signal 
can  be  used  to  obtain  temperature.  (2)  Higher  up,  for  the  range  of  altitudes  where  the  fine-structure  multiplets  of  atomic 
oxygen  are  in  LTE  and  the  electronically  excited  intermediate  state  'D2  is  relaxed  primarily  by  spontaneous  emission  of  a 
photon,  the  Stokes  and  anti-Stokes  backscattered  signal  can  be  used  to  obtain  the  atomic  oxygen  density  and  local 
temperature.  (3)  Still  higher  up,  for  the  range  of  altitudes  where  the  fine-structure  multiplets  of  atomic  oxygen  are  not  in 
LTE  and  the  electronically  excited  intermediate  state  'D2  is  relaxed  primarily  by  spontaneous  emission  of  a  photon,  the 
Stokes  and  anti-Stokes  backscattered  signal  can  be  used  to  obtain  the  density  of  the  3P2  and  3P|  multiplet  levels  of  the 
ground  electronic  state  of  atomic  oxygen.  For  a  ground-based  instrument  a  simulation  with  3  km  range  gate  is  used  to 
show  that  the  relative  error  of  temperature  measurements  from  100  to  250  km  could  be  less  than  30%.  It  is  pointed  out  that 
this  technique  has  the  potential  of  providing  unique  data  that  addresses  the  modeling  of  satellite  drag  and  the  effects  of 
space  weather  on  the  upper  atmosphere.  In  addition,  this  technique  may  also  permit  the  detection  of  the  thickness  of  the 
temperature  inversion  layers  as  well  as  their  temperature  and  density  perturbations. 

©  2005  Elsevier  Ltd.  All  rights  reserved. 

Keywords:  Middle  and  upper  atmosphere;  Atomic  oxygen  density;  Temperature;  Resonant  Raman  scattering;  Remote  sensing 


1.  Introduction 

The  upper  middle  atmosphere  and  the  lower 
thermosphere  contain  at  the  lower  boundary  the 
coldest  part  of  the  atmosphere  and  at  the  upper 
boundary  abrupt  transition  to  the  higher  tempera¬ 
tures  of  the  thermosphere.  The  steep  rise  in 
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temperature  is  accompanied  by  many  interesting 
chemical  and  dynamical  effects.  Prominent  among 
the  chemical  effects  is  the  steep  rise  in  the  atomic 
oxygen  density  and  its  dominant  role  in  cooling  the 
atmosphere  (Sharma  and  Roble,  2001,  2002). 
Prominent  among  the  dynamical  effects  is  the 
breaking  of  gravity  waves  and  energy  deposition 
in  the  mesosphere  and  lower  thermosphere  leading 
to  large  perturbations  in  temperature  and  density. 
Existing  range-resolved  remote-sensing  techniques 
for  diagnostics  of  temperature  and  density  in  the 
mesosphere  are  based  on  lidar,  operating  on  metal 
atom  transitions,  dependent  on  the  backscattered 
radiation  from  meteoric  metal  layers  (Na  (Gardner, 
1989),  K  (von  Zahn  and  Hoffner,  1996),  and  Fe 
(Gelbwachs,  1994)).  Although  these  lidar  techniques 
are  confined  to  below  120  km,  they  provide  a  wealth 
of  information  on  temperature,  winds  and  density 
perturbations.  The  technique  proposed  in  this  paper 
addresses  the  critical  gap  of  range-resolved  tem¬ 
perature  data  in  the  region  above  that  limit.  At 
those  altitudes,  there  are  a  number  of  resonance 
scattering  lidar  techniques  that  addressed  the 
measurements  of  minor  atmospheric  species, 
e.g.,  excited-state  N2  (Garner  and  Dao,  1995), 
aurorally  excited  species  (Collins  et  al.,  1997),  and 
metastable  He  (Gerrard  et  al.,  1997).  In  contrast, 
the  technique  discussed  in  this  paper  is  based  on 
resonant  Raman  scattering  and  targets  the  measure¬ 
ments  of  temperature  by  probing  a  major  atmo¬ 
spheric  component. 

We  have  recently  suggested  a  Raman  lidar 
experiment  for  in  situ  measurement  of  atomic 
oxygen  density  and  temperature  (Sharma  and 
Dao,  2005).  This  experiment  determines  the  O 
density  and  temperature  in  the  immediate  vicinity 
of  the  spacecraft  in  the  lower  thermosphere.  In  this 
article  we  extend  this  experiment  to  remote  sensing 
from  the  ground.  In  the  experiment  proposed  here, 
we  would  resonantly  excite  the  0(3P2)  ->  (3D2)  and 
0(3P|)-»(iD2)  transition  by  two  lasers  tuned  to 
630.0304  and  636.3776  nm,  respectively.  While 
resonant  Raman  lidar  has  been  studied  and 
discussed  for  lidar  measurements  of  atmospheric 
molecules  (Hochenbleicher  et  al.,  1976;  Measures, 
1984),  resonant  Raman  for  atoms  has  not  been 
generally  considered  for  lidar  measurements.  It  is 
generally  accepted  that  atoms  has  negligible  Raman 
cross  sections.  The  resonant  Raman  differential 
cross  section  of  O,  as  shown  later,  is  10-29- 
10_26cm2sr-1 — comparable  to  the  molecular 
Raman  differential  cross  section  (Measures,  1984). 


The  backscattered  Raman  signal  for  the  laser  tuned 
to  630.0304  [636.3776]  nm  is  detected  at  636.3776 
(anti-Stokes  transition)  [630.0304  (Stokes  transi¬ 
tion)]  nm,  while  the  atom  undergoes  transition 
!D2^3P,  (!D2->3P2),  giving  a  signal  proportional 
to  the  3P2  [3Pj]  density  in  a  given  space  time  interval. 
Up  to  about  300  km  altitude  the  rate  of  decay  of  the 
0(‘D)  atoms  is  dominated  by  collisions  with  N2  and 
02.  The  scattering  cross  section  up  to  these  altitudes 
is  therefore  inversely  proportional  to  the  square  of 
the  density  of  N2  and  02.  The  ratio  of  the  Stokes 
and  anti-Stokes  signal  however  depends  only  on  the 
ratio  of  the  densities  of  the  two  fine  structure  levels, 
[0(3P2)]/[0(3Pj)].  Up  to  at  least  350  km  altitude, 
these  two  fine-structure  levels  are  in  thermodynamic 
equilibrium  (Sharma  et  al.,  1994),  the  ratio  of  the 
densities  of  the  two  fine  structure  levels  gives  the 
local  temperature  T,  using  the  relation  [0(3P2)]/ 
[0(3P 1)]  =  (5/3)  exp(227.7/7)  where  227.7  K  is  the 
spacing  between  the  two  levels  (Fig.  1).  For  altitudes 
greater  than  about  300  km,  when  the  0(*D)  atoms 
decay  mostly  by  emission  of  a  photon,  and  the  fine 
structure  levels  are  still  in  thermodynamic  equili¬ 
brium,  the  Stokes  and  anti-Stokes  backscattered 
Raman  signals  give  us  the  local  temperature  as  well 
as  the  atomic  oxygen  density. 

The  technique  described  in  this  paper  relies  on  the 
moderate  differential  cross  section  for  resonant 
Raman  scattering  in  the  region  above  90  km,  where 
the  spectral  line  shape  is  dominated  by  Doppler 
broadening.  We  begin  by  calculating  the  resonant 
cross  section  for  the  two  Raman  transitions.  Our 
discussion  of  the  principles  of  temperature  measure¬ 
ment  is  followed  by  an  assessment  of  the  technique’s 
sensitivity  and  error  analysis.  As  mentioned  earlier, 


J=1 

J=2 


Fig.  1.  Energy  level  diagram  relevant  to  the  Raman  lidar.  The 
transition  from  3P2  to  ‘D2,  with  630.0304  nm  radiation,  results  in 
the  Stokes  line  at  636.3776  nm.  The  transition  from  3P2  to  'D2 
with  636.3776  nm  results  in  the  anti-Stokes  line  at  630.0304  nm. 
The  3P0  level,  98.8  K  (68.7  cm-1)  from  the  3P|  level,  is  not  shown 
because  it  is  not  involved  in  the  resonant  Raman  excitation. 
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the  temperature  is  derived  from  the  relative  signal 
strength  of  two  transitions  from  ground  states 
which  are  separated  by  an  energy  gap  of 
227.712K.  We  then  use  a  prototypical  setup  to 
project  the  performance  of  such  a  lidar  system.  The 
engineering  aspects  of  implementing  the  lidar 
experiment  are  not  dealt  with  in  this  paper. 

Hochenbleicher  et  al.  (1976)  in  a  laboratory  study 
for  the  applicability  of  resonance  Raman  scattering 
for  remote  analysis  of  air  pollutant  molecules  “found 
that  the  influence  of  the  absorption  of  the  scattering 
gas  can  reduce  the  observed  resonance  Raman  to 
normal  Raman  scattering  levels.”  These  authors  were 
dealing  with  densities  at  sea  level  involving  a  large 
number  of  collisions  in  a  radiative  lifetime.  To  obtain 
measurable  signals  they  had  to  deal  with  strong 
transitions  (short  radiative  lifetime)  so  that  a  photon 
would  be  emitted  before  the  excited  state  could  be 
collisionally  deactivated.  Strong  transition  in  turn 
induces  self  absorption,  negating  the  advantage  of 
resonance  Raman  cross  sections.  We  are  dealing  with 
about  six  orders  of  magnitude  smaller  density  at 
higher  altitudes.  With  reduced  collisional  deactiva¬ 
tion,  we  can  use  much  weaker  transitions.  We  can 
also  deal  with  moderate  transitions  under  the  thresh¬ 
old  of  self  absorption  and  take  full  advantage  of 
enhanced  resonance  Raman  scattering  cross  sections. 
At  higher  altitudes  when  the  collisional  deactivation 
of  the  excited  state  is  negligible  compared  to  the 
radiative  deactivation  the  resonant  Raman  cross 
section  becomes  even  larger.  When  the  product  of 
the  dipole  moment  operators  in  the  numerator  is 
about  the  same  as  those  involved  in  the  radiative 
decay  (Eq.  (2)),  as  is  the  case  under  consideration, 
and  the  excited  state  decays  only  by  emission,  the 
resonant  Raman  cross  section  becomes  independent 
of  the  strength  of  the  transition  to  the  intermediate 
state.  The  effects  described  here  are,  as  the  theory 
outlined  below  shows,  are  quite  general  and  not  due 
to  any  special  set  of  circumstances. 

2.  Raman  scattering  cross  section 

The  differential  cross  section  for  Raman  scatter¬ 
ing  is  given  by  (Kramers,  1964;  Mizushima,  1970; 
Loudon,  1983) 

du  =  axol  {/j(Ss-H)|f)(tt]E-»llfl 

d Q  c4  4"  \(£»  -  E,  -  hoj  +  i(hrv/2)) 

(/’I(g-P)|a)(»l£s -UMA  2  m 

(Eb  -  Ef  +  hco)  )  ’  () 


where  hco  is  the  energy  of  the  incident  photon,  hcos 
the  energy  of  the  scattered  photon,  |  i)  the  wave- 
function  of  the  initial  state  of  the  atom,  jo)  the 
wavefunction  of  the  intermediate  state  of  the  atom, 
\f)  the  wavefunction  of  the  final  state  of  the  atom,  p 
the  dipole  moment  vector  operator,  Tv  the  decay 
constants-1  (radiative  and  nonradiative)  of  the 
intermediate  atomic  state,  e  is  the  polarization 
(electric)  vector  of  the  incident  photon,  and  es  the 
polarization  (electric)  vector  of  the  scattered 
photon. 

Since  we  are  dealing  with  resonant  Raman 
scattering  for  which  Ev  =  2?,-  +  hco,  the  contribution 
of  one  intermediate  state  is  expected  to  be  much 
larger  than  that  of  any  other  intermediate  state  and 
in  addition  the  contribution  of  the  first  term  in 
Eq.  (1)  is  expected  to  be  much  larger  than  that  of  the 
second  term.  We  may  therefore  write 


d<x 


coco: 


d  Q 


{/j(ss  • »»)!«) 


hrv/2 


(2) 


Following  Zare  (1987),  we  write  the  following: 
/An\  »/2  m=+  1 

£-fl=  T  kE4W,)  (3) 

\  ^  /  Af=-1 

and 


/4ti\1/2 

■  I*  —  I  ,  j  2_,  Esm  Y im{Qp),  (4) 

where  p0  (£o)>  M-i  (£-i)>  /h  (ei)  are  the  projections  of 
the  vector  p  (e)  on  the  unit  vectors  along  the  z,  (1/ 
2)l/2(x—iy),  and  —(1/2 )1/2(x  +  iy)  axes,  respectively. 
Radiation  polarized  along  x-axis  is  obtained  by  the 
coherent  (summing  before  taking  absolute  square) 
combination  ( 1  /2)1  ^2(e_  i— e  i)  while  that  polarized 
along  y-axis  is  represented  coherently  by  the 
combination  (i'/(2)1/2)(e_i  +  ei). 

In  our  lidar  experiment  we  take  the  incident  beam 
traveling  in  the  e_i  direction  with  its  polarization 
vector  along  the  z-direction  ( M  =  0)  (Sharma  and 
Levin,  1973;  Sharma,  2004).  To  distinguish  the 
backscattered  radiation  from  the  Rayleigh  scattered 
radiation,  we  will  investigate  only  Stokes  (red 
shifted)  and  anti-Stokes  (blue  shifted)  Raman 
scattered  photons.  We  are  therefore  restricted  to 
the  case  where  the  initial  and  final  states  are 
different.  The  backscattered  light  can  be  polarized 
either  in  the  z-direction  (m  =  0)  or  the  £]  direction 
( m  =  1).  To  determine  atomic  oxygen  density  and 
temperature  we  choose  the  initial  and  final  states 
from  the  three  multiplet  levels  of  the  ground  state  of 
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atomic  oxygen  (0(3P/),  7  =  0,  1,  and  2)  and  the  only 
intermediate  state  is  0(!D2).  This  intermediate  state 
is  connected  to  all  three  7  levels  of  the  ground  state 
by  an  electric  quadrupole  transition  but  is  con¬ 
nected  only  to  7  =  2  at  630.0304  nm  and  7  =  1  at 
636.3776  nm  by  a  magnetic  dipole  transition  (http:// 
physics. nist.gov/cgi-bin/ASD/linesl. pi)  with  Ein¬ 
stein  A  coefficient  0.00563  and  0.00182  s-1,  respec¬ 
tively  (Table  1).  Since  the  electric  quadrupole 
transition  probabilities  for  these  two  transitions 
are  at  least  two  orders  of  magnitudes  smaller  than 
the  magnetic  dipole  transition  probabilities,  we  will 
ignore  the  former.  The  electric  quadrupole  transi¬ 
tion  probability  for  the  0('D2)->0(3Po)  transition 
at  639.1733  nm  is  8.60xl0-7s-1 — a  number  so 
small  that  we  will  ignore  it  and  will  not  consider  the 
transitions  from  or  to  the  O(3P0)  level  further. 

The  Einstein  A  coefficient  for  the  emission  of  a 
photon  with  a  given  direction  of  polarization  for  the 
0(’D2)->0(3P7)  is  given  by  (in  Gaussian  units) 

AC Dj  -*•  3P/)  =  ^  (^5^)l(3PyllMmagll1D2)|2, 

(5) 

where  the  expression  in  the  brackets  is  the  reduced 
magnetic  dipole  moment  matrix  element. 

The  values  of  the  reduced  dipole  matrix  elements 
are 

K3P2llAtmaEllID2)|  =  |{1D2||/rmag||3P2)| 

=  6.70  x  10-23  (6) 

and 

©1/2 

K'D2||^mag||3P,)| 

=  3.87  x  10-23.  (7) 

Using  the  Wigner-Eckart  theorem  we  can  write 

(3P/Wl(£,  •  P)l'D2,m)(lD2m|£oMo!3Py,m) 

=  (— l)Me0£_ji/C(712;  w0m)C(217;  m,  M,  rri) 

x  (3P/'||^mas||1D2)(1D2|l/imag||3P/)5m-jm+j^,  (8) 


where  m,  M,  and  m'  are  magnetic  quantum  numbers 
of  the  initial  atomic  state,  scattered  photon,  and 
final  atomic  state,  respectively.  The  Cs  are  Clebsch- 
Gordan  coefficients  and  5 m\m+M  is  a  Kronecker 
delta.  Substituting  the  relation 

C(J\2;m,  0,m)C(217';m,  M,m  +  M) 

=  [5(27'  +  1)]1/2  W(J\  17';  If) 

f 

x  C(1  If;  0,  M,  M)C(JJ'f;  m,  -M  -  m,  -M),  (9) 

where  W  is  a  Racah  coefficient,  into  Eq.  (8)  and 
inserting  the  result  into  Eq.  (2)  which  upon 
averaging  over  the  initial  magnetic  quantum  num¬ 
bers  and  summing  over  the  final  ones,  gives 

do-  _  coco]  '  *  N  2 
d Q=~?~ 


— )  |<3Pr|Lumagll‘D2} 


x(1D2||^mag||3i>/)|2  x  Yj  E  ( 

/=  1,2M=0,1  \ 

x  W2(J1 17';  2/)C2(l  If;  0,  M,  M). 


5(27'  +  1) 
27+1 


(10) 


M  denotes  the  polarization  vector  of  the  scattered 
beam  ( M  =  0  for  scattering  without  change  in  the 
direction  of  polarization  and  M  =  1  for  scattering 
with  change  in  the  direction  of  polarization), /is  the 
rank  of  the  tensor  mediating  the  scattering  process. 

The  0(*D2)  state  of  atomic  oxygen  is  rapidly 
quenched  by  02  and  N2;  the  recommended  quench¬ 
ing  coefficients,  by  the  Sander  et  al.  (2003)  are  3.2 
(+ 1.2)  x  10-n  exp(70/7)cm3s-1  and  1.8(+1.2)x 
10-n  exp(l  10/7)  cm3  s-1;  the  value  at  200  K  being 
4.5(+1.2)x  10-ucm3s-1  and  3.1(±  1.2)  x  10-n 
cm3s-1,  respectively.  Taking  the  density  of  02  and 
N2  at  100 km  altitude  to  be  1.6xl012cm-3  and 
6.5  x  1012cm-3,  respectively,  one  estimates  the 
decay  probability  T„  =  273.5  s-1.  Figs.  2  and  3  plot 
the  densities  of  the  major  thermospheric  species  and 
temperature,  respectively,  as  function  of  altitude 
(Hedin,  1991).  Because  of  its  dependence  upon 
density,  rv  decreases  with  altitude  and  therefore 
differential  Raman  scattering  increases  with  altitude 
reaching  a  limiting  value  when  rv  =  0.00745  s-1,  the 


Table  1 

Spectroscopic  parameters 


Transition 


Wavelength  (nm) 


Einstein  coefficient  (s  ’)  x  103 
magnetic  dipole  moment 


Einstein  coefficient  (s  ')  x  105 
electric  quadrupole  moment 
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Fig.  2.  Number  density  of  atomic  oxygen  in  0(3P()  (red),  0(3P2) 
(blue),  atomic  oxygen  total  (black),  molecular  oxygen  (dash 
green)  and  molecular  nitrogen  (dash  black)  as  a  function  of 
altitude.  The  density  of  fine  structure  level  0(3Po)  is  not  shown. 


Fig.  3.  Temperature  as  a  function  altitude. 


spontaneous  radiative  decay  constant.  It  is  likely 
that  the  rate  coefficient  for  quenching  of  0(1D2)  by 
ground  state  atomic  oxygen  0(3P),  not  found  in  the 
literature,  is  small  because  it  would  require  dumping 
almost  2eV  of  electronic  energy  into  translational 
motion.  We  have  therefore  chosen  to  ignore  it. 
The  differential  Raman  scattering  cross  section 
(cm2sr_1)  at  100  km  altitude  becomes 


£)  =  ^(4-81  x  1049)|  (3P /'  |  |//mag  1 1 1 D2) 

x(1D2||Mmag||3P/)|2 

X  E  E  (^T^\  W\j\\j'\2f) 


/=1,2M=0,A  2/+1 

xC2(ll/;  0 


(ID 


For  scattering  without  change  in  the  polarization 
vector  ( M  =  0)  and  7VT,  the  Racah  coefficient  and 
Clebsch-Gordan  coefficients  are  both  nonzero  for 
only  /  =  2;  hence  only  f  —  2  contributes  to  the 
scattering  process.  The  scattering,  in  other  words,  is 
mediated  by  a  tensor  of  rank  2,  and  the  expression 
for  the  scattering  cross  section  at  100  km  altitude 
becomes 

(S)100||  =  ^(4‘81  X  1  °49)  I  (3p/ 1  lA*mas  1 1 1 D2 ) 

x(,D2|ffimag||3Py)|2x 

x  W2(J\  17';  22)C2(1 12;  0, 0, 0).  (12) 

For  3P2-+3Pi  Stokes  Raman  transition  (J  =  2, 
J' =  1,  g>=  15,867.862  cm-1,  o)s  =  15,709.597 

cm-1)  without  change  in  the  direction  of  the 
polarization  vector  of  the  incident  beam,  the 
resonant  Raman  scattering  cross  section  at  100  km 
altitude  is  (d(T/d(2)100  2_>1j|  =  6.2  x  10-22  cm2  sr"1. 

For  scattering  with  change  in  the  polarization 
vector  (M  =  1)  and  J'J,  the  Racah  coefficient  and 
Clebsch-Gordan  coefficients  are  both  nonzero  for 
/  =  1  and  2;  hence  both /  =  1  and  f  —  2  contributes 
to  the  scattering  process.  The  expression  for  the 
scattering  cross  section  (Eq.  (11))  at  100  km  altitude 
becomes 


c4 


(4.81  x  1049)|{3P/'||/imag||1D2) 


x(1D2||/imag||3P.,)|2 


xC2(ll/;  0,1,1). 


W2(J\\J')  2f) 

(13) 


For  the  3P2-+3P[  Raman  transition,  substituting 
the  values  of  Clebsch-Gordan  and  Racah  coeffi¬ 
cients,  we  get 


'dcr\ 

df2/  10o,2_>.ii_l 


cocol 

c4 


(4.81  x  1049)  x  0.1 


x  |<3PiM/xmag||1D2> 
x  (1D2||//nag||3P2>|2 


(14) 


giving  (d(r/df2)100>2_^i  x  =  2.8  x  10“21  cm2  sr-1.  The 
total  cross  section  for  Raman  scattering 
630.0304  nm  photons  into  636.3776  nm  photons  at 
100km  altitude  is  therefore  (d<r/dfi)100j2^.I;i+H  = 
3.4  x  10-21  cm2  sr-1.  When  Doppler  broadening  is 
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included,  the  peak  cross  sections  are  respectively 
1.6  x  1<T27  and  1.9  x  10-27cra2sr-1. 

Proceeding  in  a  similar  manner  the  differential 
cross  section  for  the  ^  ->3P2  Raman  transition,  for 
Raman  scattering  636.3776  nm  photons  into 
630.0304  nm  photons  at  100  km  altitude,  is 


100, 1  — »2,_L-t- 1|  \u“/  100,1— *2,11  \uiV  100,  l->-2,  X 

=  1.0  x  10~21  +4.7  x  10~21 
=  5.7  x  10~21  cm2  sr*1. 


3.  Principles  of  temperature  measurement 

The  method  for  deriving  temperature  from  lidar 
Raman  signals  is  best  illustrated  with  the  use  of  Fig. 

1  which  shows  a  partial  energy  level  diagram,  not 
drawn  to  scale,  corresponding  to  the  Raman 
transitions.  The  principles  we  applied  are  similar 
to  that  demonstrated  by  Gelbwachs  on  Fe  atoms  in 
the  mesosphere  (Gelbwachs  1994).  In  the  diagram, 
the  energy  separation  (in  Kelvin  units)  between  the 
7  =  2  and  the  7  =  1  levels  is  227.712  K.  The  ratio  of 
population  at  7  =  2  level  to  7  =  1  level  in  thermo¬ 
dynamic  equilibrium  at  temperature  T  K  is  given  by 
the  Maxwell-Boltzmann  relationship 

—  =  (5/3)exp(227.712/T).  (15) 

ti\ 

With  a  temperature  of  200  K  at  100  km,  the 
Boltzmann  factor,  right-hand  side  of  Eq.  (15),  is 
estimated  at  5.2.  Because  the  population  in  the  7  = 

2  level  is  5.2  times  greater  than  the  7  =  1  level  and 
the  Stokes  (from  7  =  2  to  7  =  1).  Since  the  Stokes 
differential  cross  section  is  about  0.6  times  the  anti- 
Stokes  cross  section,  its  volume  backscatter  coeffi¬ 
cient  is  ~3  times  higher  than  that  of  the  anti-Stokes. 
Therefore  if  the  oxygen  atoms  are  in  thermal 
equilibrium  at  the  atmospheric  temperature  the 
ratio  of  Stokes  to  anti-Stokes  backscatter  signal 
can  be  used  to  measure  the  temperature. 

The  number  of  photons  recorded  by  the  lidar  is 
given  by  the  lidar  equation  as 

S  =  >i(E/hv)Az^Tno^,  (16) 

where  r]  is  the  composite  efficiency  of  the  system,  E 
is  the  transmitted  laser  pulse  energy,  hv  the  photon 
energy,  A z  the  range  gate,  A  the  receiver  area,  z  the 
altitude,  n0  the  density  of  O  atoms  and  dcro/dfi  the 
differential  cross  section  for  the  Raman  transition. 


The  laser  line  width  and  wavelength  stability  are 
assumed  to  be  less  than  the  Doppler  width  (0.8  GHz 
at  200  K).  The  efficiency  rj  includes  the  round  trip 
atmospheric  transmission,  the  transmission  through 
the  optics,  and  the  quantum  efficiency  of  the 
photomultiplier  detector.  The  lidar  equation  is  used 
to  estimate  the  photon  counts  for  the  test  cases.  In 
comparing  the  signal  for  the  anti-Stokes  wavelength 
(2^5  -  636nm:7  =  1  ->  7  =  2)  to  that  for  the 
Stokes  wavelength  (2$  =  630nm:7  =  2  7  =  1), 

one  notes  that  all  parameters  in  the  lidar  equation 
stay  constant  except  photon  energy,  number  den¬ 
sity,  differential  cross  section  and  atmospheric 
transmission.  The  ratio  of  Stokes  to  anti-Stokes 
signal  is  related  to  temperature  as  follows: 

Ss  2s  do~2i/dfi»2  / j rj\ 

S’ as  2as  dcn/dO  «i 

In  Eq.  (17),  the  subscript  S  stands  for  the  Stokes 
transition,  AS  for  anti-Stokes,  d<X2i  /dO  is  the 
differential  cross  section  for  the  Stokes  transition 
from  7  =  2  to  7  =  1  level,  and  nt  and  n2  are  the 
number  densities  of  O  atoms  in  the  7=1  and  7  =  2 
levels,  respectively.  Assuming  that  the  laser 
pulse  energy  is  independent  of  wavelength,  the 
signal  ratio  under  condition  of  thermodynamic 
equilibrium  among  the  fine-structure  levels  of 
atomic  oxygen  is 

5s  is  d»J./da  9, . /227.712  K', 

=  )■  (18) 
In  Eq.  (18),  #/is  the  degeneracy  of  level  7(=  27  +  1) 
and  227.712  K  is  the  energy  separation  in  .  Kelvin 
units  (=  AE/ks).  Using  Eq.  (11)  for  the  differential 
cross  section  and  noting  that  the  number  of 
scattered  photons  is  inversely  proportional  to  the 
cube  of  the  wavelength,  we  can  write 


—  =  -3^exp 
OAS  2as 


227.712  K 
T 


Near  the  altitude  of  100  km,  the  temperature  is 
related  to  the  ratio  of  the  Stokes  and  anti-Stokes 
signals  R  (=  Ss/Sas)  by 

„  227.712  K 

T  = - , — .  (20) 

ln((2AS/2s)3R) 

Differentiating  the  expression  of  T  in  Eq.  (20) 
relative  to  R,  we  can  show  that  the  rate  of 
temperature  change  is  as  follows: 


AT  = 


227.712K  R 
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The  sensitivity  factor,  defined  as  the  ratio  of 
fractional  change  of  T  over  fractional  change  of  R, 
is  ~1  for  the  conditions  near  the  altitude  of  100  km. 
If  we  assume  that  the  change  in  R  is  dominated  by 
the  change  in  the  J  —  1  population  or  change  in 
anti-Stokes  signal,  a  10%  temperature  detection 
limit  requires  measuring  about  8%  change  of  the 
anti-Stokes  signal.  That  corresponds  to  a  SNR  of  10 
or  an  anti-Stokes  signal  of  100  photon  counts  in  the 
photon-limited  regime  and  zero  background  noise. 
Eq.  (21)  predicts  that  the  sensitivity  'factor  will 
improve  in  the  thermosphere  because  of  higher 
temperatures  as  well  as  nearly  equal  populations  of 
the  J  =  1  and  J  =  2  fine-structure  levels. 

A  more  general  expression  of  the  temperature 
error  can  be  derived  as  follows.  If  one  can  neglect 
sky  irradiance  background  and  detector  noise,  the 
uncertainty  in  R  will  be  dominated  by  the  shot  noise 
of  photon  counts  and  the  temperature  error  is  given 
by  the  expression 

\T  T 

■7-  =  2277T2K'/(1/'Ss)  +  (I/^)'  (22) 

According  to  Eq.  (22),  error  in  T  is  inversely 
proportional  to  the  number  of  backscattered 
photons  and  hence  to  the  square  root  of  the 
integration  time. 

It  should  be  noted  again  that  the  measurement  of 
T  is  independent  of  atmospheric  density.  At  lower 
altitudes  when  the  collisional  decay  of  0(‘D)  atoms 
is  faster  than  the  spontaneous  radiative  decay,  both 
stokes  and  anti-stokes  Raman  cross  sections  are 
pressure  and  temperature  dependent.  The  tempera¬ 
ture  dependence  is  due  to  the  temperature  depen¬ 
dence  of  quenching  rate  coefficients  by  collisions 
with  02  and  N2  molecules.  Since  the  pressure  and 
temperature  dependences  cancel  out  in  the  ratio,  the 
temperature  measurement  is  insensitive  to  tempera¬ 
ture  and  density  change.  That  makes  temperature 
remote  sensing  a  robust  technique  in  the  upper 
mesosphere  and  lower  thermosphere  even  though 
both  temperature  and  atmospheric  density  vary 
considerably. 

At  higher  altitudes  when  the  decay  of  0(‘D) 
atoms  is  due  only  to  spontaneous  radiative  decay 
and  the  fine-structure  levels  of  O  are  still  in 
thermodynamic  equilibrium,  the  backscattered  sig¬ 
nal  no  longer  depends  upon  the  local  N2  and  02 
density  and  this  experiment  can  provide  both 
atomic  oxygen  density  and  temperature.  Density 
measurements  are  beyond  the  scope  of  this  paper. 


4.  Error  analysis  of  lidar  measurements 

The  signal  to  noise  ratio  of  lidar  backscatter  with 
the  wavelength  tuned  in  the  Stokes  resonance 
(As  =  630nm:7  =  2  -*■  J  =  1)  and  the  anti-Stokes 
resonance  (Aas  =  636  nm:  J  =  1  -*  J  —  2)  is  com¬ 
puted  with  the  use  of  Eq.  (16),  the  lidar  equation. 
The  assumptions  about  the  transmitter,  receiver, 
and  other  factors  affecting  the  efficiency  of  the  lidar 
are  summarized  in  Table  2.  To  calculate  the  overall 
efficiency,  the  following  values  have  been  assumed: 
one-way  atmospheric  transmission  of  0.65,  telescope 
efficiency  of  0.9,  filter  transmission  of  0.65,  and 
detection  efficiency  of  photomultiplier  detector  of 
0.20.  In  this  first-order  estimate  of  the  SNR,  errors 
in  cross  section  due  to  laser  tuning  errors  and  laser 
line  width  errors  are  considered  negligible.  In  the 
first  test  case,  the  signal  strength  for  a  ground-based 
lidar  with  characteristics  as  mentioned  is  displayed 
in  Fig.  4.  The  atmospheric  model  used  in  the 
calculations  is  the  Mass-Spectrometer-Incoherent- 
Scatter-Extended  90  (MSISE  90)  (Hedin,  1991]  with 
thermospheric  temperature  approaching  964  K. 

To  measure  temperature  in  the  lower  thermo¬ 
sphere,  the  signals  are  integrated  for  lh.  The 
estimated  relative  errors  of  the  temperature  profile 
are  shown  in  Fig.  5.  Measurements  from  a  ground- 
based  instrument  with  accuracy  better  than  30%  are 
possible  for  altitudes  between  120  and  250  km.  This 
technique  compliments  the  existing  remote-sensing 
techniques,  the  upper  limit  of  this  temperature 
technique  has  been  extended  from  the  altitudes  of 
minor  metallic  species,  ~  120  km,  to  the  lower 
thermosphere  near  300  km.  This  has  the  potential 
of  providing  unique  data  that  addresses  the  model¬ 
ing  of  satellite  drag  and  the  effects  of  space  weather 
on  the  upper  atmosphere. 


Table  2 

Lidar  parameters  used  in  the  SNR  estimates 


Parameter 

Ground  lidar 

Aperture  diameter 

3.5m 

Laser  energy  (J/pulse) 

0.05 

Repetition  rate  (Hz) 

30 

Range  gate 

3  km 

Quantum  efficiency  of  photon  counting  (PMT) 

0.2 

One  way  transmission 

0.65 

Optical  efficiency: 

Interference  filter 

0.65 

Telescope 

0.9 

Overall  efficiency  rj 

0.049 
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ALTITUDE  (km) 


Fig.  4.  Signal  strength  in  photon  counts  per  pulse  as  a  function 
of  altitude  for  a  ground  based  lidar. 


ALTITUDE  (km) 


Fig.  5.  Relative  error  in  temperature  for  a  ground-based  lidar 
with  signal  integrated  over  1  h.  Table  1  lists  the  lidar  parameters. 


5.  Summary 

We  have  calculated  the  resonant  cross  sections  for 
the  Stokes  and  anti-Stokes  atomic  oxygen  and 
introduced  a  new  lidar  technique  for  temperature 
measurement  in  the  90-290  km  region.  The  lidar 
wavelengths  are  resonant  with  the  anti-Stokes 
transition  (Aas  =  636  nm:7  =  1  ->  J  —  2)  and  the 
Stokes  transition  (As  =  630  nm:/  =  2  ->  J  =  1)  for 
atomic  oxygen.  At  100  km,  the  temperature  relative 
error  is  approximately  equal  to  the  error  of  the  ratio 
of  the  anti-Stokes  over  Stokes  signals.  Based  on  the 
prototypical  setup  we  project  that  a  lidar  can  be 
built  to  routinely  measure  range-resolved  tempera¬ 
ture  in  the  region  where  no  alternative  diagnostics 
exist.  This  system  could  be  significant  to  the  study 
of  gravity  waves,  neutral  density  drag  and  space 


weather  effects  in  the  upper  mesosphere  and  lower 
thermosphere. 
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of  altitudes  where  the  fine-structure  multiplets  of  atomic  oxygen  are  in  LTE.and  the  electronically  excited  intermediate  state  is  relaxed  primarily  by 
spontaneous  emission  of  a  photon,  the  Stokes  and  anti-Stokes  backscattered  signal  can  be  used  to  obtain  the  atomic  oxygen  density  and  local  temperature. 
(3)  Still  higher  up,  for  the  range  of  altitudes  where  the  fine-structure  multiplets  of  atomic  oxygen  are  not  in  LTE  and  the  electronically  excited  intermediate 
state  is  relaxed  primarily  by  spontaneous  emission  of  a  photon,  the  Stokes  and  anti-Stokes  backscattered  signal  can  be  used  to  obtain  the  density  of  the  3P2 
and  3Pi ,  multiplet  levels  of  the  ground  electronic  state  of  atomic  oxygen.  For  a  ground-based  instrument  a  simulation  with  3  km  range  gate  is  used  to  show 
that  the  relative  error  of  temperature  measurements  from  100  to  250  km  could  be  less  than  30%.  It  is  pointed  out  that  this  technique  has  the  potential  of 
providing  unique  data  that  addresses  the  modeling  of  satellite  drag  and  the  effects  of  space  weather  on  the  upper  atmosphere.  In  addition,  this  technique 
may  also  permit  the  detection  of  the  thickness  of  the  temperature  inversion  layers  as  well  as  their  temperature  and  density  perturbations. 
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